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Abstract 

Existence of self-similar, superposed pulse-train solutions of the nonlinear, 
coupled Maxwell- Schrodinger equations, with the frequencies controlled by 
the oscillator strengths of the transitions, is established. Some of these ex- 
citations are specific to the resonant media, with energy levels in the config- 
urations of A and iV and arise because of the interference effects of cnoidal 
waves, as evidenced from some recently discovered identities involving the 
Jacobian elliptic functions. Interestingly, these excitations also admit a dual 
interpretation as single pulse-trains, with widely different amplitudes, which 
can lead to substantially different field intensities and population densities in 
different atomic levels. 
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I. INTRODUCTION 



The generation of shape preserving excitations in nonlinear media has been the subject 
of extensive research in diverse areas of physics, ranging from hydrodynamics [1], particle 
physics [2] to quantum optics [3-5] and optical communications [6]. Starting from the ex- 
planation of the solitary waves in shallow water [7], these solutions of the nonlinear wave 
equations have found ample experimental verifications. In optical context, the nonlinear na- 
ture of the coupled Maxwell- Schrodinger equations [8], describing the interaction of classical 
radiation with matter in a resonant media, has naturally generated tremendous interest in 
the study of the pulse [3,9-13] and pulse-train solutions [14-17]. Recently the continuous 
pulse-train soliton solutions have been observed experimentally for the two-level system [18]. 
For a detailed overview, the interested readers are referred to Refs.( [19,20]), apart from the 
above references. 

The common factor that governs the existence of the self-similar excitations in various non- 
linear systems is the critical balance between nonlinearity and dispersion [21]. This severely 
restricts the set of solutions, in most of the cases, to combinations of the Jacobian elliptic 
functions [22] and also leads to definite relationships between the amplitudes, widths and 
velocities of these modes. The periodic elliptic functions, characterized by the modulus pa- 
rameter m and smoothly interpolating between hyperbolic secant and tangent to periodic 
cosine and sine functions, for m values one and zero respectively, describe localized pulses 
for m = 1 and pulse-trains for < m < 1. 

In this light, the recent findings of a number of identities involving the superposition of 
elliptic functions [23-25] and their application to various nonlinear systems [26] makes it 
extremely interesting to find the implications of the same for the multi-level resonant media. 
A number of phenomena, like self [27] and electromagnetically induced transparencies [28], 
pulse sharpening [29], pulse cloning and dragging [30], just to name a few, are ascribable to 
the existence of shape-preserving excitations in resonant media. A great deal of attention, 
both analytical and numerical, has been paid to the cases in which the resonant atoms are 



2 



asymptotically in the ground or excited state. Pulse train solutions, characteristic to excited 
media, have been studied, for two- level [14,15] and for three and five-level [16] configurations. 
Well-known techniques for generating solutions e.g., inverse scattering methods [21,31,32] 
and Backlund transformations [33] have also been employed for higher level atomic media, 
albeit under restrictive conditions on the propagation constants and other parameters of 
these systems. Further progress, in the understanding of the dynamics of the inhomoge- 
neously broadened three level system, has recently been achieved by the derivation of an 
area theorem [34] . 

In this paper, we first show the existence of novel pulse train solutions, specific to the A 
and TV-type, nonlinear, resonant media, taking recourse to the above mentioned identities, 
involving elliptic functions. Interestingly, these solutions can either be viewed as a linear 
superposition of cnoidal waves or as single cnoidal waves, with widely different, modulus 
dependent amplitudes for different pulse trains. We find that, for some of these exact 
solutions, the modulus parameter is controlled by the oscillator strengths of the atomic 
transitions, as compared to other pulse train solutions found so far, where m appears as 
a free parameter. Interference effects, originating from the superposition nature of these 
solutions, lead to substantially different field intensities and population densities in different 
atomic levels. We then point out the existence of pulse train solutions in the four-level case, 
similar to the ones found in other multi-level systems. 

II. SUPERPOSED PULSE TRAINS IN A SYSTEMS 

We start with the three-level A system, because of its wide applicability and then proceed 
to the four-level case, with the energy levels in the configuration of N, from which the other 
lower ones will follow under limiting conditions. Shape preserving solutions, in the form 
of superposition of a pair of different cnoidal waves (like a superposition of sn(x, m) and 
cn(x, m)), have been obtained earlier by Hioe and Grobe [16], for the three and five level 
systems. These exact solitary waves can have a variety of shapes, because of interference. 
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The solutions to be discussed here are superpositions of an odd number of cnoidal waves of a 
given type (no combinations involving, say both sn(x, m) and cn(x, m)), with appropriately 
displaced arguments. As will be seen later in the text, these superposed solutions can have 
widely different amplitudes. This can be understood from the fact that, both cn(x, m) and 
sn(x, m) functions take values ranging from —1 to 1, whereas the dn(x, m) function has only 
positive values. Hence, the superposed solutions involving cn(x, m) and sn(x, m) functions 
can have a much smaller amplitude as compared to the one involving dn(x, m) function, 
whose amplitude can be substantially higher. 

For the A system consisting of the ground state i, excited state e and the interme- 
diate state /, we assume the fields to have slowly varying pulse envelopes: E a (z, t) = 
£ a (z,t)e~ 1 ^ at ~ kaZ ^ + c.c. . In the rotating wave approximation, £ e , and £f act on tran- 
sitions ei and ef, respectively. The Rabi frequencies, Q a = 2^^, with d a being the dipole 
matrix element of the a-th transition, are also slowly varying functions of space and time. 
We assume, for simplicity that, all the fields are resonant with their respective transitions. 

In slowly varying envelope approximations, the resonant coupled Maxwell-Schrodinger 
equations are 

d Qf Q e 

i——C e = Ct Cj , 

dr 2 s 2 " 

%— — C* = Cp , 
dr 2 ' 

i%-C f = f -C e (1) 

dr t 2 K J 

d 

—Q f = in f C e C} , 

= tfi e c e c; . (2) 

In above, we have neglected the relaxation terms, since the pulse widths are taken to be 
smaller than the relaxation times. C a (a = e, i, f) represents the probability amplitude of 
finding the atom in the state \a >. 

The parameter fi a is given by fi a = AnAf ^j-r 2 = 27re „^" where N is the density of the 



atoms and f a being the oscillator strength for the transition of frequency u a . The atomic 
system is assumed to be continuously distributed in a non-dispersive host medium. The 
coordinate r = t — | measures time relative to the pulse center and ( = z. 

We look for self-similar solutions of Eqs. (1) and (2) i.e., solutions depending upon a 
single variable, X = (q( — Ft). Here, is the pulse duration and, as will be seen later, 
q and T will be related via the pulse velocity and another parameter, characterizing the 
solutions. The superposed pulse-train solutions, to be discussed here, appear in both A 
and N type media and can be made to satisfy a wide range of initial conditions, like no 
occupancy of the desired atomic levels to partial occupation of all the levels. It should be 
pointed out that partial occupation of levels necessarily implies initially induced coherence 
in the atomic system. We note that, the pulse train solutions of A, V and other odd-level 
systems [16], consisting of a matched pair of elliptic functions, do not satisfy the nonlinear 
equations of the A" system. 

It can be shown that, the following ansatz solutions: 

d = S , C f = b f C , C e = b e D , Q f = A f S , Vl e = A e C ; (3) 

satisfy the Maxwell-Schrodinger equations, provided the functions C, S and D are as defined 
below and the constant coefficients bf, b e and A e , Af are appropriately related. C, S and D 
are linear superpositions of the Jacobi elliptic functions cn(X, m), sn(X, m) and dn(X, m), 
respectively: S = E?=iSn[X + 4(i ' 1)K(m) , m] and C and D are analogously defined. Here, 
p is an odd integer, K(m) is the complete elliptic integral of the first kind and m is the 
modulus parameter. The above solutions are possible for odd integral values of p, since the 
cross terms in the right hand side of Eqs.(l) and (2) are cancelled because of the identities 
of the type [23], 

si(o?2 + cfa) + cyclic permutations = , 
ci(i~2 + S3) + cyclic permutations = , 

Ci(o?2 + ^3) + cyclic permutations = . (4) 
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Here, si = sn(X, m) , s 2 = sn(X + ^El, m) and s~ 3 = sn(X + m) and other functions 

are similarly defined. Although, for definiteness, henceforth we consider only the p = 3 case, 
the consistency conditions on the parameters, written below, are identical for all values of 
p. Introducing an additional parameter /i = 2qTm, for future convenience, one finds, 

\A e \ 2 = 4T 2 m^ , 

\A f \> = ^mp; ~ 1] , (5) 
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be = 



h = ~^\ (6) 

and also the constraint, /i e = jif. The positive definite character of the pulse intensity, 
requires that fj, e > fj, . 

The superposed character of the solutions, for higher values of p, leads to significant 
differences between different pulse-train and electron amplitudes, as will be explicated later. 
It also manifests in the conservation of probabilities for the electrons. For the p = 3 case, 
J2 a \C a \ 2 = 1 leads to, 



K "^- 1 + |[(^+ 1 ) 2 -«flF 
Here, q = dn(|K(m), m) and satisfies, 



(7) 



q A + 2q 6 - 2(1 - m)q - (I - m) = . (8) 

The values of q range from one to zero continuously, when m varies between zero and one. It 
can be checked that, when, m — 1, the conservation law is identically satisfied, making the 
above equation vacuous; in this case one obtains the pulse solutions of the A system. For 
m — 0, right hand side takes value one, which is ruled out on physical grounds. For other 
values of m, the above equation can be numerically solved to obtain the range of values 
of jj, for a given /i e . It should be noted that, since fi e > A*, the variable l± ^ i takes values 



between zero and one. It is found that for the superposed solutions to exist, the condition 
0.59 < < 1 should be satisfied. As depicted in Fig.l, 
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FIG. 1. Graph of the modulus parameter m versus x = 

the modulus parameter m can lie between one and zero, for the above values of Hence 
the choice of m (0 < m < 1) determines \i through Eq.(7), which in turn yields qT, through 
\i = 2qTm. Hence the solutions are characterized by the two free parameters m and Y. It 
can be easily seen that D has a much higher amplitude as compared to C and S; hence 
the |e > state has a higher population. The superposition nature of the solutions enables 
one to have these widely differing amplitudes. It can be checked that exchange of C and D 
also leads to allowed solutions. In that case the \f > state population can be made large as 
compared to \i > and |e > states. 

For the purpose of comparison, we consider the non-superposed p = 1 case. Taking 
Ci = 6jSn(X, m) and replacing D and C by corresponding cnoidal waves one finds, after 
appropriately choosing the fields, that, the parameter relationships as obtained in Eqs.(5,6) 
are unchanged. However, the probability conservation now yields substantially different 
result: 



(9) 
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It is clear from the above expression that, for, bi — 1, no pulse train solutions are possible, 
as m — 1, for this case. Since /x//x e lies between zero and one, it follows from Eq.(9) 
that, \bi\ 2 takes values ranging from m and 1. It should be pointed that, like the previous 



case, < m < 1. One is then finally left with three free parameters, T, m and bi, in their 
appropriate range of values. It is possible to obtain another solution by exchanging dn(X, m) 
with cn(X, m), in the above ansatz. This choice leads to the constraint, 

^ = (i-V ■ do) 

A*e m 

III. PULSE TRAIN SOLUTIONS OF THE N SYSTEM 

In case of the N system, one needs to consider seven coupled equations, due to the 
presence of an additional energy state \v >, in between the states |e > and \f >, of the A 
system. Here, one has the additional parameter, [i v = 2ne ^* v , f v being the corresponding 
oscillator strength. These nonlinear equations have been recently found to possess propa- 
gating pulse solutions [13]. Introducing, the additional electron amplitude C v = b v D and 
the Rabi frequency fl v = A V C in the ansatz solutions, we only write below, the new and 
modified parametric relations, as compared to the A system: 

2qTm = n v , (11) 



\A V \ 2 = AY 2 m\l 
|A e | 2 = 4r 2 m^, (12) 



b 



V 



iqA v 



bf = • (13) 

Note that here, fi v determines 2qTm, unlike the A system. The positive definite character 
of the pulse intensities now requires the inequalities, 

A*/ > Ate > fJ. v ■ (14) 
In the present case, the conservation of probability yields, 



^-i + ^[(^ + i) 2 -w] ' 1 j 

This leads to the constraint on the oscillator strengths: 0.59 < < 1, for < m < 1. It 

is worth noting that, as compared to A-system, in the present case the oscillator strengths 
determine modulus parameter m. 

The uniform group velocity, as determined by our solutions is given by, 

- = - + ^-. (16) 
v c 2T 2 m V ; 

As the width of the pulse train becomes shorter (r gets larger), the group velocity approaches 
the background medium velocity; however, for smaller values of m, the pulse train velocity 
can be significantly smaller than the medium velocity. We have checked that a different 
superposition type solution, e.g., Ci = S and exchange of C and D for other fields, do not 
alter the above conclusions. However, the population of different levels, as well as the field 
amplitudes can be substantially different, from the earlier case. It is straightforward to see, 
from Eqs.(ll) and (15), that unlike the A-system, the A^-system is characterized by one free 
parameter Y. 

It is worth noting that, for the two level case, corresponding to fx e — fj, v — fj,f, as well as 
for the V system, for which fj, e — fj, v , the constraint equation originating from the probability 
conservation yields, 

m = l+4(q 2 + q). (17) 

Since < m < 1 and q > 0, it is easy to see that, the above equation can not be satisfied. 
Hence these type of excitations are exclusive to A and A" type systems. 

A better understanding of the above solutions can be provided by a dual interpretation 
of these excitations, through a set of generalized Landen transformations, found recently. 
For the p = 3 case, the transformation formula yields, 

, / ~ \ r / \ / 4K(m) , 8K(m) . , 

adn(x, m) = [dn(ax, m) + dn(ax H — , m) + dn(ax H — , m)J , (18) 
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where a = [dn(0, m) + dn(^, m) + dn(^, m)] and 

(1-g) 2 

m — m T-, —tt; , (19) 

(l + g) 2 (l + 2g) 2 v ; 

with q = dn(|K(m), m)g. For the sn(x, m) and cn(x, hi) functions, the formulae have same 
widths in right hand side as in dn(x, hi); however, the amplitude a is replaced by f3, given 
by 

n r / x /4K(m) , ,8K(m) , , 

(3 = [cn(0, m) + cn( K J , m) + cn( K \ m)] . (20) 

o o 

The facts that, dn(x, m) takes only positive values and cn(x, m) have values ranging from 
+ 1 to —1, can be seen to be the reason behind the large value of a as compared to that of 
f3. Hence, the superposed cnoidal waves can have widely varying amplitudes, which leads to 
significant differences in atomic level occupations. 

IV. PURE CNOIDAL SOLUTIONS FOR THE N SYSTEM 

Since the four-level system has not been systematically analyzed for pulse train solutions, 
we would like to point out that, like the two-level case, pure cnoidal waves, with comparable 
electron and field amplitudes, can also be obtained as solutions to the N system. As has been 
mentioned earlier, matched pair type pulse train solutions of the odd-level atomic systems 
[16] are not solutions of the four-level system. It can be seen that, the following ansatz 
solutions: 

Ci = b iS n(X, m) , C f = b f dn(X, m) , C c = b c cn(X, m) , C v = b v dn(X, m) , (21) 

and 

Q f = A/sn(X, m) , Q e = A c dn(X, m) , fi v = A v dn(X, m) ; (22) 

satisfy the Maxwell-Schrodinger equations, provided \bi\ 2 < 1 and the constant coefficients 
bi,bf,b e ,b v and A e ,Af,A v are related. The consistency conditions on all the parameters are 
now modified: 
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2qTm 
lb~F 



|A,| 2 = 4r 2 [i-^] , 

\A f \ 2 = 4T 2 m[^ - 1] , 

fJ> v 
2 



\aj 2 = 4r 



2 t^e 



L e I rrj - ) 



be 



ifi v A e bi 
iA v bi 



b 
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2T ' 
A e * m/if ' 



[1 + L_^(^ff__^)] . (23) 



\h\ 2 m A*/ 

For a given ground state occupancy, the modulus parameter m is determined by Eq.(23). 

The velocity is now modified and depends upon \bi\ 2 : 

1 = 1 

v c 2T 2 m ' V ; 

For the purpose of considering limiting cases of the iV system, we observe that, when the 
modulus parameter m equals one, the cnoidal functions go over to hyperbolic ones, thereby 
yielding localized pulse solutions. The pulse trains are obtained, when < m < 1, the m = 
value being ruled out in the present case. Since Me ^" is positive, it is clear from Eq. (23) 
that, for the N system, |6j| 2 must be less than one, for the latter type of solutions. However, 
for the two level case, where /i e — ii v — iif, the above restriction does not apply. In that 
case, bi = 1 and one obtains the cnoidal solutions, similar to a set of solutions, obtained 
by Crisp and Arecchi et al. in Ref. [14], in an inhomogeneous media. Exchanging dn(X, m) 
with cn(X,m) in the above ansatz, leads to |6j| 2 = m; in this case one obtains the solutions, 
analogous to another set of solutions given in Ref. [14], under inhomogeneous conditions. 
This is also clear from the analysis of the A system, which is obtained from the N system 
under the condition, fj, e = fj,f, with \i = /i v . The above point also indicates that for the A 
system, obtained as a limiting case, \bi\ 2 < 1, in order for the pulse train solution to exist. 
One obtains localized pulse solutions for |6j| 2 = 1, since in that case m — 1. For the V 
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system, for which \x e = the conclusions are similar to the two-level case. A comparison 
with Ref. [16] shows that the above solutions form a subclass of the matched pair solutions, 
for the A and V system. This result is expected since matched pair type superposition 
solutions are not possible for the N system. 

For the purpose of better appreciation of the relative amplitude variations of the super- 
posed solutions, as compared to the pure cnoidal ones, Fig. 2 depicts two of the pulse trains 
with their cnoidal constituents. The widely different amplitudes of the superposed dn(X, m) 
and sn(X, m) type cnoidal waves and their respective constituents are clearly visible. 
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FIG. 2a. Plot depicting the superposed D(X,m) (solid curve) and its three 
dn(X, m) constituents (dotted curves), for m = 0.7. 
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FIG. 2b. Plot depicting the superposed S(X, m) (solid curve) and its three 
cnoidal constituents (dotted curve), for m — 0.7. 

Fig. 3 shows, the superposed solutions with their non-superposed counterparts, clearly 
bringing out their differences. 
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FIG. 3a. Plots depicting the superposed solution D(X,m) (solid curve) of Eq.(3) and its 
non-superposed counterpart (dotted curve) in the same units, for m = 0.7. These plots 
represent the electron amplitudes in the state |e > of the A system. 
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FIG. 3b. Plots depicting the superposed sn(X, m) solution (solid curve) of Eq.(3) and its 
pure cnoidal counterpart (dotted curve), for m = 0.7 and 6, = 0.9869. With b { = 1 these 
plots represent the field strengths Qf of the A system, in the same units. 

Fig. 4 reveals the same for different parameter values. One sees that, for m = 0.3 the 
amplitude of superposed sn(X,m) solution attains a very small value. It is extremely in- 
teresting to observe that, the excited state probability density can be made very small and 
other densities enhanced, through these novel solutions, which is not possible through pure 
cnoidal waves. It is clear that, this desirable feature of the pulse trains owes its origin to 
the superposition nature of the solutions. 
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FIG. 4a. Plots depicting the superposed D{X,m) (solid curve) and its non-superposed 
counterpart (dotted curve) in the same units, for m = 0.3. 
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FIG. 4b. Plots depicting the superposed sn(X, m) solution (solid curve) and its pure 
cnoidal counterpart (dotted curve), for m = 0.3 and 6j = 0.9568 . 

It is natural to enquire the effect of detunings on the above solutions and to find if the 
previous methods of incorporating detuning (particularly the case A e = A f in the A system 
[12]) would also work for deriving the pulse train solutions. 

V. CONCLUSIONS 

In summary, we have found that delicate interference phenomena involving cnoidal waves 
give rise to certain novel shape-preserving excitations, exclusive to A and A" type multi-level 
atomic media. It was also observed that matched pair-type self-similar solutions, possible 
in odd-level systems, do not occur in the even-level A" system, although a variety of simple 
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pulse-train solutions are possible in the latter case. Interestingly, the superposition nature 
of the above mentioned excitations of the nonlinear media, makes it possible for different 
atomic levels to have widely different population densities and different pulse trains possess 
significantly differing amplitudes, features not present in simple pulse-train solutions. For 
the iV-system, the frequencies of the propagating pulses, which govern the spatio-temporal 
behavior of the atomic population densities, are not arbitrary, as have been the case so far, 
with the pulse-train solutions of various systems. The modulus parameter m, which ap- 
pears in the velocities and also controls the period K{m) is determined by the propagation 
constants. In light of the potential usefulness of the multi-level systems, starting from infor- 
mation storage [35,36] to quantum computation [37], we hope that, some of these desirable 
features of the excitations may find practical applications. It is interesting to observe that, 
presence of additional non-linearities leads to chirping of the pulse train of the two-level 
systems [38]. The effect of similar non-linearities on the superposed and non-superposed 
pulse trains solutions, for the three and four-level systems is worth further investigation. 

One of the authors (PKP) acknowledges useful discussions with Prof. A. Khare. 
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